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Abstract 

A simulation study of the magnetic electron drift vortex (MEDV) mode turbulence in a magneto- 
plasma in the presence of inhomogeneities in the plasma temperature and density, as well as in the 
external magnetic field, is presented. The study shows that the influence of the magnetic field inho- 
mogeneity is to suppress streamer-like structures observed in previous simulation studies without 
background magnetic fields. The MEDV mode turbulence exhibits non- universal (non-Kolmogorov 
type) spectra for different sets of the plasma parameters. In the presence of an inhomogeneous 
magnetic field, the spectrum changes to a 7/3 power law, which is flatter than without magnetic 
field gradients. The relevance of this work to laser-produced plasmas in the laboratory is briefly 
mentioned. 

PACS numbers: 52.25.Gj, 52.35.Fp, 52.50.Jm, 98.62.En 
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I. INTRODUCTION 



The generation of magnetic field in 
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aser-produced plasmas in the laboratory 
121 ] is an fascinating and rich field of research. Possible 
mechanisms to spontaneously generate magnetic fields in a plasma include the Biermann 
battery [l^ j. which is associated with non-parallel density and temperature gradients, and 



and in the universe |!8 
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the Weibel instability 



14j | . in which the electrons have a non-isotropic temperature. The 



Weibel instability may be responsible for the generation of large-scale magnetic fields in the 
Universe, as well as in inertial fusion plasmas HQ, 1 171 ]. while the Biermann battery has 
been proposed as apossible mechanism to produce mega-Gauss magnetic fields in laser- 



produced plasmas 



1|. There are recent observations of megagauss-field topology changes 



and structure formation in laser-produced plasmas 



generation of magnetic fields in plasmas [l,|3j and effects attributed to self-generated mag- 
netic fields, such as transport of energy along surfaces jl8| and the insulation of laser-heated 
electrons from the target interior {la . flo| . inspired the investigation of magnetic surface 



121 



The observations of spontaneous 



waves in plasmas 



201 ] . This model was generalized to investigate the generation of magnetic 



fields from temperature and density gradients in the plasma [211 ] . The nonlinear properties 
of MEDV modes in a nonuniform plasma with the equilibrium density and electron tempera- 



ture gradients have been investigated analytically without 



22 



231 ] , and with [24j background 



magnetic fields. The spectral and statistical properties 25|, as well as the generation of large 
scale magnetic fields 26( by the MEDV mode turbulence, have been investigate both theo- 
retically and numerically. The generation of steep, non-Kolmogorov spectra of streamer-like 
structures by MEDV turbulence has been observed in recent simulation studies 27]. 

In this paper, we report a simulation study of the MEDV mode turbulence in the presence 
of gradients in the equilibrium electron temperature and electron density, as well as in the 
background magnetic field. We concentrate our study on the magnetic field effect and its 
influence on the spectral properties of the MEDV mode turbulence. The MEDV mode 
turbulence involve a competition between nonlinear zonal flows 28|, which we define as 
nonlinear MEDV modes with a finite scale in the direction of the equilibrium plasma density 
and temperature gradients, and streamers that have a finite scale perpendicular to the 
plasma gradients. We will here assume that the background magnetic field gradient is in the 
same direction as the plasma gradients. Our work presented here ignores collisional effects. 
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The latter are discussed in Ref. 



29[ that includes frictional and thermal forces, as well as 



gradients of the unperturbed magnetic field. 

The manuscript is organized in the following fashion. The governing nonlinear equations 
for the two-dimensional (2D) MEDV modes are presented in Sec. II. The results of computer 
simulations are displayed in Sec. III., and the observed spectral properties of the MEDV 
turbulence are discussed in Sec. IV. Finally, the results are briefly summarized and discussed 
in Sec V. 



II. NONLINEAR EQUATIONS 

We here derive the governing equations for nonlinearly interacting 2D MEDV modes 
in an inhomogeneous plasma containing equilibrium electron density, electron temperature 
and background magnetic field gradients. Following Ref. 22] , we will assume a 2D geometry 
in the xy plane, where the perturbed magnetic field is directed along the z axis so that 
B = B(x,y,t)z, where z is the unit vector along the z axis. The governing equations for 



the wave magnetic field B and electron temperature fluctuations Ti are given by 



22 



and 



|r (-B - —V 2 b) = ^-^{B, V 2 B} - — {n , T x } (1) 
at \m Ho en o J {l J 'O en o) mn o 

+ —{B, T + TJ + ( 7 - l)-^{n Q , B} = 0, (2) 

respectively. Here e is the magnitude of the electron charge, m is the electron mass, 
/io = l/c 2 eo is the magnetic permeability in vacuum, c is the speed of light, eo the elec- 
tric permittivity in vacuum, and 7 = 5/3 is the ratio between the specific heats. The 
equilibrium electron number density no and the electron temperature To are assumed to 
have a gradient along the x axis. In Eqs. flTJ and we have introduced the Poisson 
bracket notation 

r, 1 = df_dg _ d£dg_ 
ftn.dy dy dx ' 

We have neglected the scalar nonlinearity [22| in Eq. (TTjK which becomes important only on 



timescales much larger than the electron gyro period [22]]. On a longer 
of the scalar nonlinearity can be important, and numerical simulations 
causes dipolar vortices to gradually transform to monopolar vortices. 



timescale the effect 



301 ] indicate that it 
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Dividing the magnetic field as B = Bq + Bi, where Bq is the large-scale background 
magnetic field and B\ is the perturbations, and assuming that the equilibrium quantities 
n , T and B depend only on the coordinate x, we have 

^-(-B 1 -^—V 2 B 1 ) =^—({B ,V 2 B 1 } + {B 1 ,V 2 B 1 })-—{n ,T 1 } (4) 

and 

^ + —({B , T\} + {B u Ti}) + (7 - l)^{n , B{\ = 0, (5) 
ot /i en A^oeno 

where we have used that {B , T } = {no, B } = {B , V 2 B } = 0. 

In order to investigate the basic properties of the nonlinear system of equations (j3J) 
and ((Sj), it is convenient to scale it with its typical length and time scales. Noting that 
the characteristic length scale of the system is the electron skin depth A e = c/uj pe , where 
w w = (noe 2 /eom) 1 / 2 is the electron plasma frequency, and that the typical time scale is the 



plasma inhomogeneity length scale divided by the thermal speed of the electrons 22j, it is 
possible to cast the system (jlJ)-(jHJ) into the dimensionless form 

|( Bl -V^)H Bl ,V^} + KB V^-f, (6) 

and 

m &r x dB x 

— = - { B,T l} - KB —-a—. (7) 

where we have normalized the spatial coordinates x and y by the electron skin depth A e , 
the time by X e /VTe{\K>\^n) 1 ^ 2 , the magnetic field by (m/e)Vre(|«|«n) 1 ^ 2 /Ae, and the tem- 
perature fluctuations T\ by (\K\K n Y' 2 T , where Vr e = (Te/m) 1 ^ 2 is the electron thermal 
speed. The normalized background plasma gradients are given by K n = (n' /n )(c/ujp e ), 
K = [(7 ~~ l) n o/ n o — (To/T())](c/uj pe ), and the normalized magnetic field gradient by 
Kb = (B' Q /B )(c/ujp e ), where the primes denote differentiation with respect to x. The 
coordinate system is chosen such that K n > 0. With this normalization, a = +1 for k > 
and a = — 1 for k < 0. Hence, the only parameters in ([6]) and ((Tj) are a and Kb, where a 
only takes the values +1 or —1. 

Linearizing the system of Eqs. and ([70, and assuming that Bi and T\ are proportional 
to exp(ik x x + ik y y — iut), we obtain the linear dispersion relation 

(uj — K B ky)[uj + k 2 (uj — K B k y )} = k 2 a. (8) 
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where u and k = xk x + yk y are the frequency and wave vector, respectively. Equation (JHJ) 
has solutions of the form 



k,, 



to 



k b (1 + 2k 2 ) ± Jk% + 4(1 + k 2 )a 



(9) 



2(1 + A; 2 ) 

For a — —1, we have unstable MEDV modes for 4(1 + k 2 ) > k 2 b when k y ^ 0. For a = +1, 
we have only stable MEDV modes, and one can find zero-frequency waves when k = 1/kb, 
in addition to the zero frequency zonal flows at k y = 0. In the limit Kb —* we retain the 
previous result [2l|, \22 \ 



h "... 



" = ± VTW^ (10) 

which predicts stable MEDV modes for a > and purely growing MEDV modes for a < 
0. The unstable case a < corresponds to a situation where density and temperature 
gradients are in the same direction and (Tq/Tq) > (7 — l)n' /rio. This instability gives rise 
to the generation of magnetic field fluctuations and is related to the first order baroclinic 
(V^o x VTi) effect, which shows the importance of the temperature fluctuations for the 
instability to take place 21] . 

The nonlinear system possesses the conserved energy integral 



£ 




T 2 ' 

{VBf + ±1 
a 



dx dy. 



(11) 



We note that the total energy is independent of k B - For a = +1 the energy integral is positive 
definite, and thus does not allow the growing of large amplitude magnetic fluctuations from 
small-amplitude noise. For a = — 1, however, the energy integral is non-definite and allows 
large-amplitude waves to grow from the linear instability discussed above. 

III. NUMERICAL STUDY 



We have adapted the nonlinear fluid code that was developed to study the evolution 
of MEDV modes in 271] in the presence of equilibrium electron density and temperature 
gradients, to also include a background gradient in the magnetic field. Consequently, the 
energy stored in these gradients excite linear as well as nonlinear instabilities in a different 
manner, than that described in 27J, because the nonlinear mode coupling interactions are 



significantly modified by the presence of both plasma and magnetic field gradients. Our 
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numerical code employs a doubly periodic spectral discretization of magnetic field and tem- 
perature fluctuations in terms of its Fourier components, while nonlinear interactions are 
de-convoluted back and forth in real and Fourier spectral spaces. The time integration is 
performed by using the 4th-order Runge-Kutta method. A fixed time integration step is 
used. The conservation of energy given by fill I) is used to check the numerical accuracy 
and validity of our numerical code during the nonlinear evolution of the magnetic field 
and temperature fluctuations. Varying spatial resolution (from 128 2 to 512 2 ), time step 
(lO- 2 , 5 x KT 3 , 1(T 3 ), constant values of K n \V$/c 2 = 0.1 and {2/3)K n X - K T \ = 0.1 are 
used to ensure the accuracy and consistency of our nonlinear simulation results. We also 
make sure that the initial fluctuations are isotropic and do not influence any anisotropy 
during the evolution. Nonlinear interactions can, however, lead to anisotropic turbulent 
cascades by migrating spectral energy in either k y ~ or k x ~ modes. The preference of 
energy transfer in either of the modes is determined primarily by the background gradients 
in magnetic or temperature and density fields that, in nonlinearly saturated state, govern 
the mode coupling interactions. In our present simulations, we largely focus on two cases, 
viz (i) small amplitude and (ii) large amplitude evolution of turbulent fluctuations. These 
two cases are characteristically different from each other in terms of dominant nonlinear 
interactions. To compare the effects of the amplitude we need to compare the nonlinear 
terms (the vector nonlinearities) with the linear terms, most importantly the coupling terms 
between the magnetic field and temperature fluctuations. Hence, we can compare the non- 
linear term {Bi, V 2 B\} to the coupling term dTi/dy in the dimensionless Eq. (6). In the 
simulations (Figs. 1 and 2) we observed that the amplitudes of the dimensionless magnetic 
field and temperature fluctuations are approximately the same for both the small and large 
amplitude cases when the turbulence is fully developed. Hence we may use T\ ~ B\ and 
compare the nonlinear term {Bi,V 2 B 1 } to dBi/dy. As an order of magnitude estimate, 
the spatial derivatives are of the order 1/L, where L is the typical length scale of the fluc- 
tuations. Hence we have the order of mag nitude estimates \{B 1 ,V 2 B 1 }\ ~ (l/L) 4 ^! 2 and 
\dBi/dy\ ~ (l/L)\Bi\. One can conclude that if B\ -C L 3 then the linear coupling term 
dominates and we have a case of weak turbulence, while if B\ « L 3 or larger, then we have a 
fully nonlinear case with strong turbulence. In our simulations, the fluctuations have typical 
length-scales of L ~ 1 (corresponding to the electron skin depth in dimensional units) so 
that \Bi\ <C 1 constitutes a weakly nonlinear case while B ~ 1 is a strongly nonlinear case. 



6 




2 4 6 2 4 6 

x X 




2 x..4 6 2 x 4 6 

FIG. 1: (Color online) Evolution of random turbulent fluctuations initialized with small amplitude. 
Nonlinear turbulent interactions lead to the formation of relatively large-scale flow directed along 
the mean magnetic field in our 2D simulations. The saturated structures in B, T, V 2 B and 
B — V 2 B are shown in the figure. The numerical resolution is 256 2 , box dimension is 27T x 27T, and 
the parameters used are a = + 1 and Kb = 0.5. 

In dimensional units, B x should be replaced by {eBi/m)\ e /\YT e {\^\^n) 1 ^ 2 ] and L by L/\ e so 
that e\Bi\/m should be compared to (Vr e /A e )(|«;|«; n ) 1 / 2 (L/A e ) 3 in dimensional units. Since 
L/\ e m 1, we then have the weakly nonlinear case e\B\\jm -C (Vr e / A e )(|K|K n ) 1//2 and the 
strongly nonlinear case e\B\\jm w (Vr e /A e )(|K|K Tl ) 1//2 in dimensional units. 

The initial spectral distribution in the magnetic and temperature fluctuations comprises 
a uniform isotropic and random amplitude associated with the Fourier modes confined to a 
smaller band of wave number {k < 0.1 k max ). While spectral amplitude of the fluctuations 
is random for each Fourier coefficient; it follows a k~ x or k~ 2 scaling. Note again that our 
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FIG. 2: (Color online) When the amplitude of initial fluctuations in magnetic and temperature 
fields is large enough, characteristic nonlinear interaction modify the flows observed in Fig (1). 
The diamagnetic-like nonlinear interactions seem to dominantly suppress the flow and lead to 
small scale isotropic turbulent fluctuations. This is shown in B, T, V 2 B and B — V 2 B. 

final results do not depend on the choice of the initial spectral distribution. The spectral 
distribution set up in this manner initializes random scale turbulent fluctuations. Since 
there is no external driving mechanism considered in our simulations, turbulence evolves 
freely under the influence of self-consistent nonlinear interactions. Note however that driven 
turbulence in the context of the MEDV mode will not change inertial range spectrum to 
be described here. The driving mechanism helps sustain turbulent interactions without 
modifying the inertial range turbulent cascades. The initial isotropic fluctuations in magnetic 
and temperature fields are evolved through nonlinear fluid Eqs. (EJ) & (J7J). The dominant 
nonlinear interactions in the inhomogeneous MEDV modes are governed by zxVB-W 2 B in 
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magnetic field equation. This nonlinearity is similar to the polarization drift nonlinearity z x 
V0 ■ V V 2 0, <p being the electrostatic potential fluctuations, in a two dimensional Hasegawa 



Mima-Wakatani (HMW) model describing drift waves in inhomogeneous plasmas 3l|, |32, 



33 



341 ] . This nonlinearity characterizes Reynolds stress forces that plays a critical role in the 



formation of zonal-flows. Analogously, one can expect generation of nonlinearly generated 
flows in underlying MEDV model here. The temperature evolution, on the other hand, is 
governed by z x VT ■ V-B nonlinearity that is identical to a diamagnetic nonlinear term 
in HMW model. The role of this nonlinearity has traditionally been identified as a source 
of suppressing the intensity of nonlinear flows in drift wave turbulence. Nevertheless, the 
presence of the linear inhomogenous background in both equations can modify the nonlinear 
mode coupling interactions in a subtle manner. Our objective is to understand the latter 
in the context of nonlinear interactions mediated by the inhomogeneous B and T fields in 
MEDV modes. 

We find from our small amplitude nonlinear evolution case (see Fig 1) that the nonlinear 
interactions in the inhomogeneous MEDV modes are typically led by the Reynolds-stress- 
like nonlinearity, i.e. z x VB ■ VV 2 5, in the magnetic field equation. Consequently, the 



-li 



te flows, as 
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34|. This 



mode structures in the saturated turbulent state are dominated by the zona 
typically observed in the Hasegawa-Mima-Wakatani (HMW) model 3l|, [ 
is demonstrated in Fig. (1) where large scale structures (zonal flow-like, k y « 0) along 
the background magnetic and temperature field gradients are developed in B, T, V 2 B and 
B — V 2 B fluctuations. Thus, small amplitude nonlinear evolution, as shown in Fig. (1), 
is consistent with the HMW model of the inhomogeneous drift-wave turbulence. Our sim- 
ulations described in Fig. (1) can be contrasted with our previous work 27| on the small 
amplitude case that was studied in the absence of a background magnetic field gradient. 
In Ref 27| . we have shown that mode coupling interactions during the nonlinear stage of 



evolution leads to the formation of streamer-like structures in the magnetic field fluctua- 
tions associated with k y « 0, k x ^ 0. These structures were similar to the zonal flows 
but contained a rapid k x variations, thus the corresponding frequency is relatively large. 
The temperature fluctuations in Fig (1) of 27|, on the other hand, depict an admixture 
of isotropically localized turbulent eddies and a few stretched along the direction of the 
background inhomogeneity. Clearly, the presence of a background magnetic field gradient 
reduces the k x variation in the steady-state flow. This leads to a considerable modification 
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FIG. 3: (Color online) Evolution of anisotropic mode structures as described by k x and k y mode 
averaged over the entire turbulent spectrum in inhomogeneous MEDV turbulence. Initially, k x = 
k y . Progressive development of anisotropy k x > k y is ascribed to the presence of background 
gradients in magnetic and temperature fields for which the anisotropy angle O = tan -1 (k x / 'ky) 
deviates continually from 45°. 

in the final structures, which now appear to look like zonal-flows (and not the streamers that 
were observed in the absence of a background magnetic field). Interestingly, the large ampli- 
tude case in our simulations unravels a completely different scenario where nonlinear mode 
coupling interactions in MEDV modes are observed to suppress the steady-state flow. This 
is shown in Fig. (2). The suppression of the zonal flow-like structures can be understood in 
the context of the diamagnetic-like nonlinearity that seems to modify the nonlinear mode 
coupling interactions in the large amplitude evolution. This nonlinearity, corresponding to 
a {T, B} term in the temperature equation, becomes gradually strong enough to nullify the 
emergence of the large scale zonal-flow- like structures. Hence, the steady state structures 
in Fig. (2) appear to possess more small scale fully developed turbulent fluctuations. We 
confirm that this state is not entirely isotropic and the background gradients nonlinearly 
maintain the anisotropic cascades in inhomogeneous MEDV modes. We elucidate this point 
in the following section. 
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IV. ANISOTROPIC MEDV CASCADES 



We quantify the degree of anisotropy mediated by the presence of large scale gradients in 
the magnetic and temperature fields in the nonlinear 2D inhomogeneous MEDV turbulence. 
In 2D turbulence, the anisotropy in the k x — k y plane is associated with the preferential 
transfer of spectral energy that empowers either of the k x and k y modes. The anisotropy in 
the initial isotropic turbulent spectrum is triggered essentially by the background anisotropic 
gradients that nonlinearly migrate the spectral energy in a particular direction. To measure 
the degree of anisotropic cascades, we employ the following diagnostics to monitor the evo- 
lution of k x mode in time. The k x mode is determined by averaging over the entire turbulent 
spectrum that is weighted by k x . 



Here Q represents any of B, T, V 2 B and B — V 2 B. Similarly, the evolution of k y mode is 
determined by the following relation. 



We can define an angle of anisotropy such that = tan -1 (k x / k y ) . It is clear from these 
expressions that the k x and k y modes exhibit isotropy when k x ~ k y for which 6 ~ 45°. 
Any deviation from this equality leads to a spectral anisotropy. We follow the evolution 
of k x and k y modes in our simulations for increasing amplitude of the initial fluctuations. 
Our simulation results describing the evolution of k x and k y modes are shown in Fig. 3. 
It is evident from Fig. 3 that the initial isotropic modes k x ~ k y gradually evolve towards 
an highly anisotropic state in that spectral transfer preferentially occurs in the k x mode, 
while the same is suppressed in k y mode. Consequently, the spectral transfer in k x mode 
dominates the evolution and the mode structures show elongated structures along the y- 
direction. With the increasing amplitude of the initial fluctuations, there exists increasing 
degree of anisotropy as illustrated in Fig 3. The larger the amplitude is, stronger are the 
nonlinear interactions. Correspondingly, there exists increasing degree of disparity in the k x 
and k y modes owing primarily to the increasing depletion in k y modes as shown in Fig 3. 

The increasing angle of anisotropy (in Fig 3), with the increasing magnitude of the initial 
fluctuations, is ascribed to the generation of highly asymmetric flows in our simulations. We 
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FIG. 4: (Color online) Energy associated with the anisotropic flows that corresponds essentially to 
the nonlinear k y = modes. Figure depicts the self-consistently generated nonlinear flows lead to 
turbulent anisotropy. The growth rate of the generation of anisotropic flow is directly proportional 
to the amplitude of the initial fluctuations. 

exemplify this point by means of Fig 4 that describes evolution of the energy associated with 
the anisotropic flows led predominantly by the k y = mode. Notably, this mode is generated 
explicitly by the nonlinear interactions. It is clear from this figure that growth rate of the 
generation of the anisotropic (dominated by the k y = mode) flow is directly proportional to 
the magnitude of the initial fluctuations. Thus the growth rate is higher when the amplitude 
of the magnetic and temperature field is large, i.e. B max ~ T max = 1.0. The smaller initial 
amplitude of the fluctuations correspond to relatively weak nonlinear interactions for which 
the generation of the k y = mode is insignificant. The nonlinear interactions do not 
introduce anisotropy and hence lead to nearly isotropic turbulent cascades. This is observed 
clearly in Fig 4 (see the curve corresponding to B max ~ T max = 0.1) which is consistent 
with the corresponding curve in Fig 3. In Fig 4, each curve is normalized with it's own 
initial value to rescale all the curves on a single plot. This further enables us to make a 
vis a vis comparison between the three different cases shown in Fig 4. Also evident from 
the figure 4 is the large amplitude fluctuations that lead to the higher growth rate of the 
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FIG. 5: (Color online) Anisotropic inhomogeneous MEDV turbulent fluctuations exhibit a /c~ 7 / 3 
like spectrum. The spectrum corresponds to an intermittent state where turbulence coexists with 
anisotropic structures. 

anisotropic flows (see dashed lines whose slope increases with the increasing amplitude of 
B max and T max ). The increasing slope, i.e. growth, associated with the anisotropic flows in 
our simulations is further consistent with the increasing angle of anisotropy that is observed 
in Fig 3. 

It is noteworthy from our simulations that the stronger nonlinear interactions pile up an 
increasing amount of turbulent energy in the k y = mode. Consequently, the turbulent 
correlation length scales tend to decrease across the flows. Hence the large amplitude simu- 
lations show decorrelated flow structures in Fig 2. This physically means that the nonlinear 
interactions led by the polarization and diamagnetic like terms in the presence of background 
gradients quench the flow (that was observed along the x-direction in Fig. 1) to introduce 
reduced turbulent decorrelated structures. 

While there exists a disparity in the spectral transfer of energy corresponding to the k x and 
k y modes, the 2D volume averaged turbulent spectrum follows a k~ 7 ^ power 
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aw, as shown 



33. 



34). The 



in Fig. (5). This spectrum is steeper than that of the HMW turbulence [3Xl . 
steepness of the observed spectrum can be ascribed to the coexistence of partially anisotropic 
flows and turbulent fluctuations in the steady state MEDV mode turbulence. 
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V. SUMMARY AND CONCLUSIONS 



In summary, we have investigated the properties of the MEDV mode turbulence nonuni- 
form magnetoplasma containing gradients in the electron temperature, the electron number 
density, and the external magnetic field. We find that the influence of the magnetic field 
gradient is to suppress the streamer-like structures observed in simulations without the 
magnetic field gradient [27|]. In addition, the steep spectra changes to spectra with a 7/3 
power law in the presence of the magnetic field gradient. We also discussed the conditions 
for an instability through a first order baroclinic effect, which would lead to the generation 
of magnetic fields under conditions where the background density and temperature gradients 
are in the same direction. It is noteworthy that the anisotropic terms in our simulations 
become important when the linear terms, corresponding to the gradients in magnetic and 
temperature fields, compete with the nonlinear terms. In such case, efficient migration of 
energy takes place between the gradients and turbulent modes that primarily lead to the 
nonlinear anisotropic flows. Our study could have relevance for laser produced plasmas in 
the laboratory, spontaneously generated magnetic fields, and effects associated with self- 
generated magnetic fields such as transport of energy along surfaces and flux limitation of 
electrons, have been observed for many years. 
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